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Direct photon production in central Pb+Pb collisions at CERN-SPS en- 
ergy is calculated within the relativistic microscopic transport model UrQMD, 
and within distinctly different versions of relativistic hydrodynamics. We 
find that in UrQMD the local momentum distributions of the secondaries 
are strongly elongated along the beam axis initially. Therefore, the pre- 
equilibrium contribution dominates the photon spectrum at transverse mo- 
menta above ~ 1.5 GeV. The hydrodynamics prediction of a strong correla- 
tion between the temperature and radial expansion velocities on the one hand 
and the slope of the transverse momentum distribution of direct photons on 
the other hand thus is not recovered in UrQMD. The rapidity distribution 
of direct photons in UrQMD reveals that the initial conditions for the longi- 
tudinal expansion of the photon source (the meson "fluid") resemble rather 
boostinvariance than Landau-like flow. 



1 



The radiation of real and virtual photons has frequently been proposed as a diagnostic 
tool for the hot and dense matter created in (ultra-) relativistic heavy-ion collisions ll|J^. 
The mean free path of photons with high transverse momentum exceeds the expected source 
sizes by one order of magnitude . Hard photons thus can give insight into the early stage 
of these reactions. 

Transverse momentum dependent upper limits for direct photon production in central 
S+Au collisions at 200^4 GeV have been published by the WA80 collaboration [Q. These 
data initiated theoretical studies which showed that direct photon production is strongly 
overestimated if the photon source is thermalized (having an initial energy density of 2 — 
5 GeV/fm^) and if one assumes that it is composed of light mesons only (say vr, 77, p, u). Due 
to the low specific entropy of these particles a reasonable final-state multiplicity of secondaries 
and a maximum temperature that is consistent with the WA80 data {T^ax < 300 MeV) can 
not be achieved simultaneously (at least if the expansion is approximately isentropic). 

Most of the calculations quoted (except those of ref. however, assumed that the 
photon source is in local thermal equilibrium and that ideal hydrodynamics can be applied to 
determine the space-time evolution of the temperature. Here, we perform a calculation within 
the microscopic transport model UrQMD, which includes the pre-equilibrium contributions 
to direct photon production. This is of particular relevance in view of the fact that local 
thermalization (i.e. locally isotropic momentum distributions) in (ultra-) relativistic heavy- 
ion collisions is probably not achieved within the first few fm/c where the high-fc^ 
photons are produced. Since neither thermal nor chemical equilibrium is assumed, the 
effects of finite viscosity (i.e. finite mean free paths) |Tl| and inhomogeneous or fluctuating 
energy density distributions [l^ are included. Comparisons of the transverse momentum and 
rapidity distributions of direct photons with those calculated within various fluid-dynamical 
models are made. 



The calculations presented here are based on the UrQMD model Il3|, which spans the 
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entire presently available range of energies from SIS to SPS. It's collision term contains 50 
different baryon species (including nucleon, delta and hyperon resonances with masses up 
to 2.2 GeV) and 25 different meson species (including strange meson resonances), which are 
supplemented by their corresponding antiparticles and all isospin-projected states. 




FIG. 1. Total cross sections of vr^vr and various other meson- meson reactions as calculated 



within the UrQMD model, data from ref. [14 



The model is based on the covariant propagation of all hadrons on classical trajectories, 
excitation of resonances and strings and their subsequent decay resp. fragmentation. UrQMD 
accounts for secondary interactions: the annihilation of produced mesons on baryons can lead 
to the formation of s-channel resonances or strings. Free cross-sections for hadron-hadron 
scattering in the hot and dense nuclear matter are employed. Comparisons of UrQMD 



calculations to experimental hadron yields from SIS to SPS are documented elsewhere []T^ 
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We briefly discuss the two most important reactions for tlie creation of direct pliotons. 
Below 2 GeV center of mass energy intermediate resonance states are excited. The total 
cross section for these reactions are given by 

21 R + 1 



R=pj2,. 
TT 



(2/Mi + l)(2/Af, + l) 



X 



J- R-^MiAh^ tot 
PcM {Mr - v^)2 + % 



which depends on the total decay width Ftot, the partial decay width Tr^MiM2 and on the 
cm. energy y/s. At higher energies these processes become less important. One then enters 
the region of t-channel scattering of the hadrons. 

UrQMD predicts a rich and non-trivial resonance structure in the meson-meson scattering 
cross section. In fig. |l| (left) the total tt+tt^ cross section as calculated within UrQMD is 
compared to experimental data [|14|. Below 1 GeV the spectrum is dominated by the p 
resonance, while at higher energy the /2(1270) gets important. 
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FIG. 2. Number of vrvr collisions as a function of the rapidity difference in the three spacial 
directions (Pb+Pb, ELab = 160^ GeV, 6 = fm). 



Figure S shows the distribution of vrvr collisions with a given relative rapidity in a given 
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spacial direction (upper indices refer to the particle number): 

One observes that most coUisions occur at \Ay\ < 3. In this range the colhsion-spectrum is 
nearly isotropic, corresponding to a temperature of T = 180 ± 30 MeV. In the longitudinal 
direction, however, there is an additional component above \ Ayz\ = 3. This is due to pions 
created in the fragmentation of a color flux tube close to projectile and target rapidities. In 
p(200y4 GeV)+p collisions ^ 15% of all pions have IvcmI > 2 [1^. Note that not all of these 
collisions are forbidden by the finite formation time of secondaries since the fragmentation 
region mesons contain one of the initial valence quarks (i.e. the sources of the color field). 
Thus, even within the formation time, these mesons may interact with half of the free meson 
cross section in the UrQMD model. Obviously, these meson-meson collisions with large 
longitudinal rapidity difference are not taken into account in hydrodynamical calculations. 
However, as will be discussed below, it is just their contribution which dominates the direct 
photon spectrum above /cr ~ 1-5 GeV. 

Let us now get to the reference equilibrium calculations of direct photon production 
within hydrodynamics. Here we assume either a three-dimensional expansion with cylindrical 



symmetry and longitudinal boost invariance [16| or a one-dimensional expansion with Landau 



initial conditions (i.e. f = at time t = 0). In both cases the dynamical creation of the fluid 
of secondaries in space-time is ignored and only the expansion stage is treated. 

In the three-fluid model fluids one and two correspond to projectile and target, while 
fluid three represents the newly produced particles around midrapidity. Fluids one and two 
are coupled via source terms leading to energy and momentum exchange. These interactions 
are due to binary collisions of the nucleons in the respective fluids and are derived from 
nucleon-nucleon data. A detailed discussion of this model as well as results (pion rapidity 
and transverse momentum spectra at CERN-SPS, baryon stopping and directed flow at AGS 



and SPS, the width of the compression shock waves etc.) can be found in refs. ||8|, p!7|jl8 



In the hydrodynamical calculations the fluid of produced particles is assumed to be the 
hottest, and thus the dominant source of high-fc^- photons. It's equation of state is that of 
an ideal gas of massive vr, 77, p, and uj mesons below Tc = 160 MeV. For T > Tc "we assume 
an ideal QGP (massless, noninteracting u, d quarks and gluons) described within the MIT 
bag-model [|19|. The bag constant is chosen such that the pressures of the two phases match 
at T = Tc (-B^/^ = 235 MeV) thus leading to a first order phase transition. 

The number of emitted direct photons in each of the three phases is parametrized as 

-^^^^■^-^'^'■'(^-) ■ 

E is the photon energy in the local rest frame. In our calculations we fix as = 5'^/47r = 0.4. 
Eq. (|^) accounts for pion annihilation (tttt p'j), and Compton-like scattering (vrp nj, 
TTf] TT'y) off a p or ?7 meson (in lowest order perturbation theory). In the QGP phase, 
alternatively, quark-antiquark annihilation {qq — > (77) and Compton-like scattering of a 
quark or antiquark off a gluon {q,q + g 9, ^ + 7) are considered. Contributions from the 
Ai meson [^], as well as the effect of hadronic formfactors 0], are neglected since they are 
of the same magnitude as higher order corrections to eq. (|^), which we have also not taken 
into account. Also, the number of processes that contribute to direct photon production 
enhances the rate only linearly, while the temperature distribution enters exponentially. 
Our main interest here is not the absolute photon yield but rather the effective temperature 
of the direct photons, i.e. the inverse slope of their transverse momentum spectrum. 

The photon spectrum emitted in a heavy-ion collision is obtained here by an (incoherent) 



integration over space-time: 



d^E—-— . (4) 



d^hrpdy J d'^xd^k 
To allow for a comparison with the hydrodynamical calculations we have considered the 
same processes also in the UrQMD model, namely n^n^ p°7, vr'^'Tr'^ P^li Tt'^^p^ — > ''^^1, 
TT^p^ TT^'-f, 7r°p^ — > 7r^7, TT^T] TT^'j . Of coursc, in UrQMD these processes are considered 



explicitly (we employ the differential cross sections given in ref. 0) and are not folded 
with thermal distribution functions, nor with a thermal A/s-distribution of meson-meson 
collisions. We find that the vrp 717 and — » 717 processes are dominant in the range 
1 GeV< kx < S GeV. This is due to the fact that, in addition to the kinetic energy, also the 
mass of the p or u can be converted into photon energy. In contrast to the thermal rate @ , 
in the UrQMD model rrip and are smeared out according to a Breit-Wigner distribution. 
For the processes vrvr p7 we have, however, assumed that the p-meson in the final state 
is produced with the peak mass of 770 MeV. Since free current- quarks and gluons have not 
been included in the present version of the UrQMD model, the processes qq —>■ (77 etc. are 
not considered there. 
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FIG. 3. Temperature in the central region as a function of time in various hydro dynamical 
models, Pb+Pb-collisions (6 = fm) at CERN-SPS. 



The temperature of the newly produced particles in the central region is depicted in fig. 
T7[]. In case of a Bjorken expansion we have assumed an initial temperature of 300 MeV 



and an initial time of tq = 0.22 fm/c, while for the Landau expansion we take an initial 
longitudinal extension of 2L = 2i?pb/7cM ~ 1-4 fm. The maximum temperature in the 



three-fluid model is determined by the coupling terms between fluids one and two (and, of 
course, by the equation of state of fluid three). The produced particles cool most swiftly 
in scaling hydrodynamics, while in the case of a Landau expansion the temperature at 
midrapidity is constant until the rarefaction waves reach the center. In this latter case 
one finds the longest-living mixed phase. This is due to the fact that we assumed purely 
one-dimensional (longitudinal) expansion. 
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FIG. 4. Transverse momentum spectrum of direct photons at midrapidity (Pb+Pb, 
^Lab = 160^ GeV, 6 = fm). The lines refer to hydrodynamical calculations including a first order 
phase transition to a QGP at Tq = 160 MeV (three-fluid model (dashed), scaling expansion (solid), 
and Landau expansion (dotted)). The result of the UrQMD calculation with all meson-meson colli- 
sions (full circles) is compared to the calculation with only "thermal" meson-meson collisions (open 
circles). Crosses show the photons from w-meson decays, which are as many as those from "thermal" 
meson-meson collisions. Above ~ L5 GeV the pre-equilibrium contribution dominates. 



The various transverse momentum distributions of direct photons, cf. fig. ^, reflect the 
distinct cooling laws. The faster the cooling of the photon source, the steeper the slope of 
the photon spectrum. An exponential fit of the spectrum in the region 2 GeV< A;t < 3 GeV 
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yields "T" = 260 MeV both for the three- fluid model and for the Landau expansion, and 
"T" = 210 MeV for the Bjorken expansion. 

If we take into account only the "thermal" meson-meson collisions with rapidity gap 
lAy^l < 3, we find a similar photon spectrum in the UrQMD model as in scaling hydro- 
dynamics. The fact that the hadronic mass spectrum in UrQMD has no upper limit leads 
to a collision spectrum of the light mesons that is equally "cool" as in the hydrodynamic 
calculations, if there a first order phase transition at a critical temperature around 160 MeV 
is assumed (cf. also ref. 0). In UrQMD, a significant part of the energy density can be 
stored into heavy meson- and baryon resonances. This keeps the pressure and temperature 
low P,pT|]. We have analyzed the temperature as a function of energy density (at normal 
nuclear matter baryon density) in UrQMD in ref. [IT3 . 



One also observes that the "pre-equilibrium" meson- meson collisions with jAy^l > 3 
dominate the emission of direct photons with large transverse momenta, > 1-5 GeV. 
Obviously, this contribution has nothing to do with high meson transverse momenta ( "tem- 
peratures") or fast collective radial expansion. It is thus not included in the hydrodynamic 
calculations. Hence, the strong correlation between the (Doppler-shifted) temperature and 
the slope of the photon transverse momentum distribution predicted by the hydrodynamical 



calculations (see also refs. |p|,[7|,P,|22| ) is not seen in UrQMD. Meson-baryon and baryon-baryon 



processes may further enhance the pre-equilibrium contribution [^]. Below = 1.5 the 
photons are dominantly produced by "thermal" meson-meson collisions and — 717 decays. 
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FIG. 5. Rapidity distribution of direct photons, calculated within the various hydro dynamical 
models (for kr = 2 GeV) and UrQMD (for kr = 2.15 GeV). 



The rapidity distribution of the secondaries (or, alternatively, of the temperature) at 
early times is very different in the various hydrodynamical models |T7|J2^ . While T{r]) (rj 
denotes the fluid rapidity) is strongly peaked around midrapidity in the three-fluid and 
Landau models, scaling flow Vz = z/t and energy-momentum conservation, d^T^'^ = 0, 
imply that the pressure p(T,^b) is independent of the fluid rapidity, drjP = 0, and depends 



only on proper time, r = \Jt'^ — . For net baryon-free matter d^^T = follows. However, 
even in the first case matter is accelerated by a (longitudinal) rarefaction wave, leading to 
a broadening of the rapidity distribution. Thus, in any case, at the late freeze-out stage 
the rapidity distribution of secondaries (around midrapidity) will be more or less flat. On 
the other hand, the rapidity distribution of direct photons with transverse momenta much 
larger than the maximum temperature of the photon source (e.g. /cr ~ 2 GeV at SPS) 
offers the opportunity to measure the rapidity distribution of the hot photon source at early 
times [^. This is demonstrated in fig. ^. In the three-fluid model and in the Landau- 
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expansion the photon rapidity distribution is strongly peaked around midrapidity. It is not 
proportional to the (squared) rapidity distribution of the pions, which are emitted at the 
much later freeze-out stage. The distribution obtained within UrQMD only resembles that for 
a boostinvariant expansion, if only "thermal" meson-meson collisions are taken into account. 
The full calculation predicts a maximum of the photon-dN / dy at midrapidity. This is so 
because the photon spectrum emitted in elementary meson-meson reactions with large Ay is 
peaked around midrapidity. Nevertheless, the photon distribution is considerably flater than 
in the three-fluid or Landau models, which show a concave curvature in this logarithmic plot. 
This is a consequence of the longitudinal velocity profile ||25] induced by the (longitudinal) 
rarefaction wave accelerating the fluid that is initially at rest. The UrQMD calculation does 
not exhibit this characteristic feature. 

In summary, we have calculated direct photon production in central Pb+Pb collisions 
at CERN-SPS energy within a microscopic transport model (UrQMD). Within this model, 
secondary hadron production is described by color flux tube fragmentation and resonance 
decay. This leads to non-thermal momentum distributions, i.e. they are strongly elongated 
along the beam axis initially. The contribution from this pre-equilibrium stage to direct 
photon production is found to dominate at transverse momenta above ^1.5 GeV. Thus, 
within the UrQMD model, the slope of the photon transverse momentum spectrum is not 
directly related to the Doppler-shifted "transverse temperature" of the source, as it is in 
hydrodynamics. 

If only "thermal" meson-meson collisions are taken into account, both the transverse 
momentum and rapidity spectra of direct photons in UrQMD resemble those calculated 
within boostinvariant hydrodynamics, if for the latter one assumes an initial temperature 
Ti = 300 MeV and initial time = 0.22 fm. This is so although QGP formation is not 
assumed in UrQMD. Hence, heavy resonances and color flux tubes allow for a rather soft 
meson-meson collision spectrum even at high entropies. 
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